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Abstract: We study a discretization of N = 2 super Yang-Mills theory which possesses 
a single exact supersymmetry at non-zero lattice spacing. This supersymmetry arises 
after a reformulation of the theory in terms of twisted fields. In this paper we derive the 
action of the other twisted supersymmetries on the component fields and study, using 
Monte Carlo simulation, a series of corresponding Ward identities. Our results for SU(2) 
and SU(3) support a restoration of these additional supersymmetries without fine tuning 
in the infinite volume continuum limit. Additionally we present evidence supporting a 
restoration of (twisted) rotational invariance in the same limit. Finally we have examined 
the distribution of scalar field eigenvalues and find evidence for power law tails extending 
out to large eigenvalue. We argue that these tails indicate that the classical moduli space 
does not survive in the quantum theory. 
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1. Introduction 

Supersymmetric theories have long been of interest to particle physicists both from a phe- 
nomenological and theoretical perspective. However for many years the study of such 
systems on the lattice was problematic. Supersymmetry is typically broken at the classical 
level in such theories and this makes it difficult if not impossible to construct supersymmet- 
ric continuum limits for such theories - see, for example the reviews ||] and references 
therein. 

Recently we have developed a formalism for discretizing supersymmetric theories in 
D dimensions with Af = p x 2 D ,p = 1, 2, . . . supercharges in which one or more super- 
symmetries are preserved exactly at non-zero lattice spacing. 1 The discretization proceeds 
from a reformulation of the continuum theory in terms of twisted fields. Such a procedure 
naturally exposes a nilpotent supercharge Q and generically leads to an action which can 
be written in Q-exact form. These properties allow us to construct supersymmetric lattice 
actions which preserve the corresponding supersymmetry Several lattice models have 
been proposed starting from these twisted continuum formulations []|, ^, [7J. In this paper 

1 Strictly the construction yields a lattice theory with complexified lattice fields but theoretical arguments 
and numerical work have provided evidence that the supersymmetry is retained at the quantum level when 
the fields are restricted to be real. 
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we will investigate the discretization developed in []|] and §. In this case the resulting 
lattice theories yield an alternative geometrical reformulation of the orbifold constructions 
of Kaplan et al. |lC|j . The explicit connection between the the twisted and orbifold lat- 
tice constructions was shown recently by Unsal [|TI]] . The two constructions share the same 
naive continuum limit and preserve the same number of supersymmetries. Indeed, we show 
later that linear combinations of the scalar and rank 2 tensor supercharges occurring in 
the twisted formulation of M = 2 SYM yield, in the naive continuum limit, the conserved 
nilpotent charge of the equivalent orbifold model. Thus the two approaches to lattice su- 
persymmetry, while differing in their details, depend on essentially the same underlying 
mechanisms for preserving supersymmetry at finite lattice spacing. 

The fermion fields that appear in the twisted formulations can be naturally embedded 
as components of one or more Kahler-Dirac fields. This geometrical interpretation of the 
fermions makes it then almost trivial to discretize the theory in a way which avoids the 
notorious problem of fermion doubling. Furthermore, using the discretization prescription 
developed in ||, the model may then be studied using numerical simulation techniques 

Hi- 

In this paper we extend previous work by deriving the action of the remaining super- 
symmetries on the twisted fields and use Monte Carlo simulation to examine a series of 
Ward identities which yield information on the restoration of full supersymmetry in the 
continuum limit. The next two sections discuss the continuum twisted theory and its tran- 
scription to the lattice. We then derive the remaining twisted supersymmetries and go on 
to describe our numerical results concerning both the Ward identities and the question of 
the restoration of rotational invariance. The bulk of our simulations are conducted in the 
phase quenched approximation in which any phase of the Pfaffian resulting from integrating 
the fermions is neglected. We show numerical results indicating that this approximation 
is likely adequate in the infinite volume continuum limit. We have also conducted simula- 
tions employing temporal antiperiodic boundary conditions for the fermions which further 
suppress the phase fluctuations. 

We also show some new results on the distribution of scalar fields in this theory. We 
argue that the form of this distribution is consistent with a lifting of the classical moduli 
space via quantum fluctuations. The final section addresses our conclusions. 

2. Twisted J\f = 2 SYM in two dimensions 

As discussed in || the action of the the two dimensional N = 2 super Yang-Mills theory 
can be written in the twisted or Q-exact form 

S = /3QTr j d 2 x Qr#, 0] + 2*12^12 + Xi2#i2 + (2.1) 

where Q is a scalar supercharge obtained by twisting the original Majorana supercharges 
of the theory. The twist consists of decomposing all fields under the action of a twisted 
rotation group. The latter is obtained as the diagonal subgroup of the direct product 
of the original (Euclidean) Lorentz symmetry with an SO(2) subgroup of the theory's 
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R-symmetry. In practice this means that we should treat the Lorentz and flavor indices 
carried by all fields (and supercharges) as equivalent resulting in a 2 x 2 matrix structure 
for the fields. The twisted fields of the theory correspond to the expansion coefficients 
when this matrix is decomposed on a basis of products of gamma matrices. Explicitly for 
the supercharge matrix q this decomposition reads 

q = QI + QpTlii + <3i27i72 (2.2) 

while a similar expression for the fermions allows the fermionic content of the theory to be 
re-expressed in terms of the set antisymmetric tensor fields (77/2, ipn,, Xi2)- The action of 
the scalar supercharge on the component fields is 2 





= Vv 




= 


Qx\2 


= B\2 


QB 12 


= [<P,Xl2 




= 2^ 




2 




= \\m 






Q6 


= 



(2.3) 

Notice that supersymmetry requires the introduction of Q-superpartners with the same 



tensor structure as the fermions. Carrying out the Q- variation on eqn. 2.1 and subsequently 
integrating over the field B12 leads to the action 

S = (3Tr f d 2 x Q[<M 2 - - Ff 2 - D^Drf 

- X12&X12]- 2xi2 (£)i^2-£>20i)-^^r ? + ^[0,^]) (2.4) 

The bosonic sector of this action is precisely the usual Yang-Mills action while the fermionic 
sector constitutes a Kahler-Dirac representation of the usual spinorial action [|l4|]. This can 
seen explicitly by constructing a Dirac spinor out these twisted fields in the following way 

A ( 2*' ~ lX A (2.5) 



iip- 



■2 



It is straightforward to see that the kinetic terms in 2.4 can then be rewritten in the Dirac 
form 

X^j.DX (2.6) 
where the gamma matrices are taken in the Euclidean chiral representation 





71 = L n 72 = ■ n ( 2 - 7 ) 



2 All fields are to be thought of as antihermitian matrices in the gauge group. 



-3- 



In the same way the Yukawa interactions with the scalar fields can be written 



where 75 in this representation is 

*-(J-°i) 

Thus the on-shell twisted action is nothing more than the usual J\f = 2 SYM action in two 
dimensions. Indeed, in flat space the twisting process can be thought of as simply a change 
of variables and is hence fully equivalent to the usual formulation. 3 

This rewriting of the theory in terms of antisymmetric tensors has two primary advan- 
tages - it allows us to formulate the theory on a curved space and as we will show in the 
next section gives a natural starting point for discretization. 

Finally it is worth pointing out that the twisted theory possesses an additional U(\) 
symmetry inherited from the remaining R-symmetry of the model which is given by 





- e iQ W 






2 — ► 


e~ ia Xi2 
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-» e~ ia V 








e 2ia <j) 


6- 


-> e~ 2ia 6 



(2.10) 

This symmetry can also be preserved under discretization and ensures the absence of 
additive mass renormalizations in the theory. 



3. Lattice Action and Symmetries 

This theory may be discretized by mapping the continuum rank p antisymmetric tensor 
fields to lattice fields living on p-cubes in a hypercubic lattice and replacing derivatives 
with appropriate difference operators. 

The choice of difference operator is very important to avoid fermion doubling - one 
must replace — > if the derivative occurs in a curl-like operation and d — > D~ if 



the derivative belongs to a divergence IjTSf . Here, Dt and D~ refer to the usual forward 
and backward difference operators respectively. Here gauge covariant versions of these 
difference operators must be used. We have employed the following definitions [16]: 

D+f(x) = U^x)f(x + M ) - f(x)U„(x) 

D+f v (x) = U„(x)U(x + /i) - U(x)U^x + u) (3.1) 



^Notice that to make this correspondence and obtain a bounded Euclidean action it is necessary to think 
of <f> as a complex matrix which is (minus) the Hermitian conjugate of <f>. 
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Notice that these definitions imply that the gauge transformations of lattice scalar, vector 
and (antisymmetric) tensor fields are 

f(x) -> G(x)f(x)G\x) 
U(x) - G(x)f lx {x)G^(x + l j) 

U»(x) - G{x)f lxv {x)G\x + n + v) (3.2) 

where G(x) = e^ x ' is a lattice gauge transformation and all fields are taken in the adjoint 
representation. From a lattice perspective these transformations are very natural and 
correspond to thinking of each field as living on a link running from the origin out to a 
vertex on the unit hypercube. The definitions of the backward difference operator follow 
by taking the adjoint 

D ^Uu(x) = fvw(x)Ul(x + Ulf^x - n) 

D'Uix) = U(x)U}(x) - U£(x - fif^x - /i) (3.3) 

Notice that this lattice divergence has the merit of demoting the rank of the lattice field as 
in the continuum and is consistent with the lattice gauge transformation rules listed above. 

The final step to constructing the lattice theory is to promote each (anti)hermitian 
continuum field to a complex lattice field. This allows / and /' to transform differently 
under gauge transformations. This in turn is required if we are to construct gauge invariant 
objects on the lattice. This doubling is actually somewhat natural in a lattice theory with 
p-form lattice fields since the underlying p-cube with p > has two orientations. The 
complexification of the vector potential A^(x) has the additional benefit of allowing the 
fields U(x) and U\ x) to vary independently under the twisted supersymmetry. In the end 
we will require the final path integral be taken along a contour where U'U = I and the 
imaginary parts of the gauge field and the fermion fields vanish. This reality condition will 
allow contact to be made with the usual twisted continuum theory [||, 13|. 

The Q-exact lattice action now takes the form 

S L = PQTr^ (\vHx)[^Mx)]+x{2(^i2(x)+Xi2(x)T 12 (x^ 

x ^ 

+ \xU*)Bi2(x) + \xi2{x)B\ 2 (x) + ^l(x)D+$(x) + \%{x){D+4(x)y \ (3.4) 

This expression will also be Q-invariant if we can generalize the continuum twisted super- 
symmetry transformations in such a way that we preserve the property Q 2 = 5q. The 
following transformations do the job 

QUfj, = Vv 

= ~D+<f> 
Q(j) = 
QX12 = B u 
QB 12 = [0,xi 2 ] (12) 

Q<t> = v 

Qv = [<M (3.5) 
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where the superscript notation indicates a shifted commutator 



[<f>, Xnv] iH = 4(x)Xp>(x) - Xnu(x)<f>(x +ii + u) (3.6) 

These arise naturally when we consider the infinitesimal form of the gauge transformation 
property of the plaquette field. Notice that gauge invariance also dictates that we must use 
the covariant forward difference operator on the right-hand side of the variation. 
The lattice field strength is given by ^^{x) = D+U v {x). 

Carrying out this lattice Q-variation leads to the following expression for the lattice 
action 

S L = /3Tr^ (-Mx)Mx)} 2 - V(xMx), !?(*)] - xl 2 (s) Xi2(*)] (12) + b\ 2 (x)B 12 (x) 



+ b\ 2 (x)F 12 (x) + B 12 {x)F l2 { X y + -(D+<f>{x)?Dffix) + -D+<t>(x)(D+^(x)^ 

- xl 2 (x)B>iMx) + x\ 2 (x)DtMx)) ~ 4> ] 2 {x)D^ X i2(x) + ip\(x) D 2 X n(x) 

~ ^l(x)D+ V (x) - lvHx)D-M^+^M^M^)} M ) (3-7) 

Finally we must integrate out the multiplier fields B\ 2 and b\ 2 resulting in the term 

^Tr^^ 12 (x)t^ 12 (x) (3.8) 

This can be written 

/3Tr ( 2I -Up- Uj,) + /3Tr (M 12 + M 21 - 21) (3.9) 



where 

Up = TV (Ui{x)U 2 {x + V)U\{x + 2)u\{x)] (3.10) 



resembles the usual Wilson plaquette operator and 

M 12 (x) = U^U^U^x + l)U\{x + 1) (3.11) 

Notice that the second term vanishes when the gauge field is restricted to be unitary which 
is equivalent to requiring Iva.A^{x) = 0. In this case the action is nothing more than the 
usual Wilson gauge action. 

Having constructed this complexified lattice theory possessing the exact scalar super- 
symmetry we will subsequently truncate it to the real line by setting the imaginary parts of 
all fields bar the scalars to zero (the scalars are required to be (anti)hermitian conjugates 
of each other as in the continuum). In |l3|] we argued that such a truncation should be 
valid for the Ward identities associated with the scalar supersymmetry at least for suffi- 
ciently large (3 as a consequence of the Q-exact nature of the action. The numerical results 
presented in [13| and in this paper are consistent with this. 



-6- 



4. Twisted Supersymmetries 



It is straightforward to construct the twisted supersymmetry transformations of the com- 
ponent fields. It follows from the matrix nature of the twist that the continuum fermion 
kinetic term can be written in the form 

S F = / d 2 xTr^j.D^f (4.1) 



where corresponds to the matrix form of the Kahler-Dirac field 

* = |^ + VV7m + Xi27i72 (4-2) 

This term is clearly invariant under —* *&T l , i = 1 ... 4 where T l are the basis of products 
of gamma matrices introduced earlier in eqn. |2.2j . 

Consider first the case V = 7172- In terms of the component fields the transformation 
^ — > ^r 4 effects a duality map 



1 
2 



-X12 



X12 



1 
2 

VV ~ e ^u (4.3) 



Such an operation clearly leaves the Yukawa terms invariant and trivially all bosonic terms. 
It is thus a symmetry of the continuum action. By combining such a transformation with 
the original action of the scalar supercharge one derives a additional supersymmetry of the 
theory - that corresponding to the twisted supercharge Qi2- Explicitly this supersymmetry 
will transform the component fields of the continuum theory in the following way 





= -e^v 


Q12VV 


= —t/ivDvi 


Q12X12 




= -\{4>A] 


Q12-B12 






= "2X12 


Ql2 2 


= B12 


Ql2<t> 


= 



(4.4) 

From the Q and Q12 transformations it is straightforward to verify the following algebra 
holds 

{Q,Q} = {Qi2,Ql 2 } = ^ 
{Q,Q 12 } = (4.5) 

where 8a, denotes an infinitessimal gauge transformation with parameter c/>. This allows 
us to construct strictly nilpotent symmetries Q± = Q ± 1Q12 in the continuum theory 
corresponding to using the (anti) self-dual components of the original Kahler-Dirac field. 
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In the same way we can try to build an additional supersymmetry by combining the 
invariance of the fermion kinetic term under ^ — > ^fF 1 with the existing scalar supersym- 
metry. This effects the following transformation of fermion fields: 

f -^ 

Xl2 ~1p2 

^ 2 

^2 -> -X12 (4.6) 

However, the Yukawas and bosonic terms are only invariant under such a transformation 
if we simultaneously make the transformation <j> — > —(f). The resultant explicit action of Qi 
on the component fields is given by 



QlA! = 


V 
2 


QiA 2 = 


~X12 


Qiipi = 




Ql^2 = 


—B12 


Qixn = 


-D 2 ~4> 


Q\B\ 2 = 




Qi$ = 





9.2 = 


D\4> 


Qi<f> = 


-2^1 



(4.7) 

Similarly the supersymmetry associated with T 2 is given by 



Q 2 ^i 


= X12 


Q2A2 


_ 7] 
2 


Q2^l 


= B12 






Q21P2 


= -\[^1>] 


Q2X12 


= D{4> 


Q2B12 


= -foW 


Q20 


= 


<h\ 


= D 2 


Q2(t> 


= "2^2 



(4.8) 

Again, we can verify the following algebra holds 

{QuQi} = {Q2,Q2} = 5_- 4> 

{Q 1 ,Q 2 } = (4.9) 
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with 5_-^ a corresponding gauge transformation with parameter —(f). This allows us to 
construct yet another pair of nilpotent supercharges in the continuum theory Q± = Q\ ± 
iQi- 

It is interesting to check also the anticommutators of these new charges Q± and Q ± . 
It is a straightforward exercise to verify the following algebra holds on-shell 

{Q+,Q-} = {Q~,Q+} = o 
{Q+,Q+} = 4(Di + »D 2 ) 

{Q-,Q-} = 4,(D 1 -iD 2 ) (4.10) 

As an example consider {Q+,Q + }il>i 

{Q+,Q + } = {Q, Qx} - {Qn, Q2} + i({Qn, Qi} + {Q, Q2} (4.11) 

Using the component transformations listed above the relevant anticommutators are 

{Q, Qi}ipi = 2Ihih 
{Q,Q 2 }V>i = 2D 1 ip 2 + 2&,Xi2] 
{Qi2,Qi}V>i = 2D 2 ^i 

{Q12, Q 2 }V>i = 2£> 2 ^ 2 + [(/>, V] (4-12) 

Thus we find 

{Q+,Q+} = 21>i^i - 2D 2 ^ 2 - [</>, V] + i(2D^ 2 + 2D 2 ^i + 2[<f>, X12}) (4.13) 
Using the equations of motion 

-2D 1 ^ 1 -2D 2 ^2-[<P,V] = 

-2D 2 ^ + 2D^ 2 + 2[<f>, X n] = (4.14) 



we can easily verify the second line of eqn. 4.10 . Notice that from these new charges Q±, 
Q± we can build spinorial supercharges of the form 



0+ 
Q- 



(4.15) 



in which case the algebra given in eqn. 4.10| represents the usual supersymmetry algebra 
in a chiral basis (up to a gauge transformation). 

Finally we will be interested in Ward identities which can be derived for a general 
operator O and take the form 

<Q i O>=0 a = l... 4 (4.16) 

for each of the four supersymmetries. To get nontrivial results the operator O must be 
gauge invariant and have U(l) charge —1. 
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5. Numerical Results 



We have implemented the RHMC dynamical fermion algorithm to simulate this lattice 
theory - for details we refer the reader to fl3|| . We have examined lattices with L 2 geom- 
etry with L in the range 3 — 8 and a variety of lattice couplings. Simulations have been 
done employing both periodic boundary conditions and thermal boundary conditions corre- 
sponding to enforcing antiperiodicity on the fermions in the temporal direction. Following 
fi~3[ | we have worked in the phase quenched ensemble - the Pfaffian arising from integra- 
tion over the twisted fermions is replaced by its absolute value within the Monte Carlo 
simulation. We do, however, monitor the phase and in the next section will discuss the the 
results of reweighting our measured observables with those phase fluctuations. Typically, 
our results derive from O(10 3-4 ) trajectories for both SU(2) and SU(3) theories where a 
single trajectory corresponds to r = 1 units of classical dynamics time. 

5.1 Scalar supersymmetry 

Consider first a series of Ward identities associated with the exact lattice supersymmetry 
Q. The simplest of these corresponds to the statement < S >= — dX Q^ = reflecting the 
Q-exact nature of the twisted action. If we integrate out the twisted fermions and the 
auxiliary field B\2 we find the following expression for the partition function 

Z = pW a N./2p-N a N./2 J D ^ DUe -PSB(^,U) det ^M(^)) (5.1) 

where N s is the number of sites and Nq the number of generators of the gauge group. 
The first prefactor arises from the fermion integration while the second derives from the 
gaussian integration over the auxiliary field. From this we find the following condition on 
the mean bosonic action as a consequence of the scalar supersymmetry 

< j3S B >= 1 (5.2) 



3N G N S 

To examine this quantity in the continuum limit we must know how to scale the lattice 
coupling j3 with the number of lattice points. Clearly the physics of the model is determined 
by the dimensionless coupling fi = g 2 A where A is the physical area. If we simply equate 
this to the corresponding lattice quantity we find the scaling 

T 2 

(3 = - 

The continuum limit is thus gotten from this equation taking L — > oo while holding \i 
fixed. Notice that the infinite volume limit corresponds then to taking the subsequent 
limit n — > oo. 

Figure [l] shows a plot of this quantity for the SU(3) gauge group for \x = 0.25 and 
/x = 5.0 as the continuum limit is approached with L = 3,4,5. While small deviations 
of order 1.5% are seen from the theoretical result based on exact supersymmetry on the 
smallest lattice with = 0.25 these appear to diminish as the continuum limit is taken. 
Furthermore, at fj, = 5.0 the data appears consistent with exact supersymmetry for all 
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SU(2) 


5*7(3) 


L 


B 


F 


B 


F 


3 


0.031(2) 


-0.032(2) 


0.13(2) 


-0.13(2) 


4 


0.0087(8) 


-0.0094(6) 


0.053(9) 


-0.055(9) 


5 


0.0032(5) 


-0.0038(6) 


0.0091(9) 


-0.0109(8) 



Table 1: QOi vs L for SU{2), SU(3) and /z = 0.25 





5?7(2) 


SU(3) 




B 


F 


B 


F 


3 


0.0334(7) 


-0.0363(5) 


0.077(3) 


-0.083(4) 


4 


0.0215(3) 


-0.0222(2) 


0.049(3) 


-0.051(2) 


5 


0.016(1) 


-0.0146(1) 


0.036(1) 


-0.0382(2) 



Table 2: Q0 2 vs L for SU{2), S77(3) and fx = 0.25 





5f7(2) 


SU(S) 


L 




F 


B 


F 


3 


0.085(3) 


-0.0808(3) 


0.21(1) 


-0.2172(5) 


4 


0.046(2) 


-0.0460(3) 


0.11(1) 


-0.1232(4) 


5 


0.029(3) 


-0.0295(1) 


0.081(6) 


-0.0790(4) 



Table 3: Q0 3 vs L for SU(2), SU(3) and fi = 0.25 



lattice sizes (the maximum deviation being 0.5%). We can easily derive other Ward iden- 
tities associated with the Q-symmetry by taking the Q-variation of the following operators 
0\(x) = rf(x)[4>(x),r](x)], 2 (x) = Xi 2 ( a; )-P 1 i2(^) and 3 (x) = ifjj l (x)D IJi 4>(x). In practice 
we have examined the integrated quantities Yl x Oi(x) and estimated the latter by selecting 
a source point x randomly on each Monte Carlo configuration. Such a procedure minimizes 
the statistical errors for a fixed amount of computation. After Q-variation we find 

Q0 2 = F^(x)F^(x) - xU D t^] 

QO s = - ^D+n - V+lVv 0] (5.3) 

The results for the expectation values of these Q-variations as a function of lattice size 
L = 3, 4, 5 are shown in tables 1,2 and 3. for \i = 0.25 and gauge groups SU(2) and SU (3). 
Similar data at fj, = 10.0 (for 577(2)) and // = 5.0 (for 577(3)) are shown in tables 4,5, 
and 6. We denote the bosonic contribution to the Ward identity by B and the real part 
of the fermionic contribution by F - the imaginary part is always small and statistically 
consistent with zero. To give a more graphical illustration of the presence of the exact 
scalar supersymmetry Q we show in figure || plots of B + F versus L for \i = 0.25, 5.0 and 
gauge group 577(3). In general we see that these Ward identities are rather well satisfied 



- 11 - 





SU{2) 


SU{3) 


L 


B 


F 


B 


F 


3 


2.7(2) 


-2.8(3) 


3.8(1) 


-3.9(1) 


4 


0.81(3) 


-0.91(3) 


1.12(4) 


-1.16(3) 


5 


0.39(1) 


-0.44(1) 


0.50(2) 


-0.55(2) 



Table 4: QOi vs L for 5^(2) Qu = 10.0) and 51/(3) (/x = 5.0) 





SU(2) 


SU(3) 


L 


B 


F 


B 


F 


3 


1.05(8) 


-1.07(8) 


1.27(4) 


-1.31(4) 


4 


0.755(9) 


-0.777(6) 


1.00(2) 


-1.003(7) 


5 


0.51(3) 


-0.53(2) 


0.67(1) 


-0.686(5) 



Table 5: Q0 2 vs L for SU(2) (fi = 10.0) and SU{3) (/i = 5.0) 





SU(2) 


SU(S) 


L 


B 


F 


B 


F 


3 


3.40(8) 


-3.24(2) 


4.51(5) 


-4.36(1) 


4 


1.87(2) 


-1.812(3) 


2.46(4) 


-2.454(2) 


5 


1.165(9) 


-1.162(2) 


1.570(30) 


-1.590(10) 



Table 6: Q0 3 vs L for SU{2) {p, = 10.0) and SU(3) (// = 5.0) 



for all lattice sizes and certainly as we move toward the continuum limit. 
5.2 The broken supersymmetries 

In this section we examine Ward identities corresponding to the supersymmetries which are 
broken after discretization. Referring to the previous section we see that these supercharges 
in general will transform a field living on one link to a neighboring link in a manner similar 
to the link constructions used in Q. This in turn implies that the variation of any closed, 
gauge invariant loop will vary into an open gauge variant loop. The latter will automatically 
have vanishing expectation value due to gauge invariance. Thus the set of operators whose 
variation under one of these link supersymmetries yields a non-trivial Ward identity is 
rather small and is further constrained by the need to obtain singlets under the additional 
U(l) symmetry described in eqn. 2.10] . We have examined the following ones 



Q12O4 = Q12 (xi 2 [<M) = \[^} 2 -xU^Xu] (5.4) 
Qi0 5 = Qi (Vl[<M) = \[^} 2 - (5.5) 
Q 2 Q 6 = Q 2 (V 2 [<M) = \[<P,W - (5.6) 
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SU(2) 


SU{3) 


L 


B 


F 


B 


F 


3 


0.077(6) 


-0.0053(5) 


0.032(4) 


-0.028(4) 


4 


0.0022(2) 


-0.0013(1) 


0.013(3) 


-0.012(3) 


5 


0.0008(1) 


-0.0004(1) 


0.0023(2) 


-0.0018(2) 



Table 7: Q 12 O a vs L for SU(2), SU{3) for /z = 0.25 





S*7(2) 


SU(3) 


L 




F 


B 


F 


3 


0.077(6) 


-0.066(6) 


0.032(4) 


-0.030(4) 


4 


0.0022(2) 


-0.018(1) 


0.013(3) 


-0.013(2) 


5 


0.0008(1) 


-0.0007(1) 


0.0023(2) 


-0.0023(2) 



Table 8: Qi0 5 vs L for SU{2), SU(3) for fi = 0.25 





SC/(2) 


SU{3) 


£ 




F 


B 


F 


3 


0.67(8) 


-0.6(1) 


0.95(3) 


-0.92(4) 


4 


0.204(7) 


-0.161(6) 


0.28(1) 


-0.25(1) 


5 


0.099(3) 


-0.072(2) 


0.124(6) 


-0.114(5) 



Table 9: Q 12 4 vs L for SU(2) (fi = 10.0) and SU(3) (fi = 5.0) 



They correspond to terms in the lattice action. Notice that they imply that the expectation 
values of different Yukawa interactions should be equal. Tables 7 and 8 show data for SU (2) 
and SU(3) and lattice sizes L = 3,4, 5 at fj, = 0.25. Similar data for SU(2) at fj, = 10.0 
and SU(3) at = 5.0 are shown in tables 9 and 10. The results for the Ward identity 
Q2OQ are very similar to those for Q\ and we omit them. Some of this data is also shown 
graphically in figures |3| and ||| which again plot B + F against L for SU{3) at fi = 0.25, 5.0 
for the operators O4 and O5. 

We see no statistically significant evidence for breaking of these Ward identities in these 
tables except for the case of Q12O4 and gauge group SU (2). In this case the breaking appear 
worse for small fi and appear to survive the continuum limit L — > 00. We conjecture that 
the explanation for these breakings lies in the use of the phase quenched approximation. In 
the next section we examine this issue more carefully and find evidence that the problem 
of phase fluctuations gets worse at small /U. Moreover, since the latter are driven by near 
zero modes of the fermion operator it is also plausible that the effect is enhanced for small 
N since the number of fermionic zero modes scales as N while the number of non zero 
modes varies as N 2 . 
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SU{2) 


SU(3) 


L 


B 


F 


B 


F 


3 


0.67(8) 


-0.64(7) 


0.95(3) 


-0.95(3) 


4 


0.204(7) 


-0.193(6) 


0.28(1) 


-0.26(9) 


5 


0.099(3) 


-0.092(3) 


0.124(6) 


-0.123(5) 



Table 10: Q1O5 vs L for SU{2) (fi = 10.0) and SU(3) (fi = 5.0) 



5.3 Phase quenching 

In this section we try to quantify the magnitt 
due to phase quenching. Initially we focus on 
tions on all fields. Tables 11 and 12 show the 
theory at fi = 0.25 and \x = 10.0 respectively. 



L 


< cos a > 


3 


0.14(1) 


4 


0.13(2) 


5 


0.14(4) 



Table 11: SU(2) fi = 0.25 



ide of any corrections to the Ward identities 
lattices employing periodic boundary condi- 
cosine of the phase < cos a > for the SU(2) 
Notice that discretization effects appear to 



L 


< cos a > 


3 


0.33(6) 


4 


0.37(1) 


5 


0.27(1) 



Table 12: SU{2) fi = 10.0 



be small and that the mean value increases with increasing fi. A plot of the distribution 
P(cos a) for SU (2) from simulations at L = 3 and fi = 0.25 is shown in figure ||. A similar 
picture for fi = 10.0 is plotted in figure ^. From this figure we see that the distribution 
possesses two peaks - one at cos a = 1 and another at cos a = —1. Thus the measured 
Pfaffian is predominately real but has an indefinite sign. From the observed weak depen- 
dence of < cos a > on L and direct observation of the distribution P(cos a) at increasing 
L we conclude that these phase fluctuations survive the continuum limit for any fixed fi. 
Notice though that the height of the cos a = — 1 peak decreases relative to the cos a = 1 
peak as fi increases which is responsible for the observed increase of < cos a > with fi. We 
conjecture that the phase becomes concentrated at cos a = 1 in the infinite volume limit 
corresponding to fi — > oo. In this limit then the phase quenched approximation would be 
exact. 

Furthermore, we conjecture that neglect of these phase fluctuation is the origin of the 
breaking of the Q12 Ward identity observed in the previous section for the SU(2) theory 
at small fi. We have attempted to check this by examining the reweighted values of the 
bosonic action. The data is shown in table 13 which shows < Sb > and its reweighted 
cousin Sg as a function of fi from simulations of the SU (2) theory with L = 3. 

While the measured value of the reweighted action has very large errors there is an 
indication that it lies closer to the exact value relative to the unreweighted action at small 
fi. Notice again though that even the naive mean action appears to approach the expected 
theoretical value for sufficiently large fi. 
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To summarize we have observed non-trivial fluc- 
tuations in the phase of the Pfaffian which appear to 
survive the continuum limit for any finite continuum 
coupling fx. These phase fluctuations appear to become 
less important in the infinite volume limit correspond- 
ing to [A — > CO. 

We conjecture that this non-trivial phase is asso- Table 13: Sb and reweighted S§ vs 
ciated primarily with the fluctuations of lowest lying A* f° r SU (2) and L = 3 
eigenvalues of the fermion operator. We have directly 

observed such states which correspond to the superpartners of the bosonic zero modes as- 
sociated with the classical moduli space. One simple way to try to reduce the problem is 
to consider the theory at finite temperature by employing antiperiodic temporal boundary 
conditions on the fermions. We consider this in detail in the next section. 

5.4 Thermal boundary conditions 

In tables 14 and 15 we show the measured values of < cos a > versus L for couplings 
fi = 0.25 and /i = 10.0 from thermal simulations using gauge group SU{2). Notice that the 
use of thermal boundary conditions does indeed push the distribution toward cos a = 1. 
This is consistent with the effect being driven primarily by the near zero modes whose 
eigenvalues are lifted to 0(1/ L) at finite temperature. The distribution for [i = 10.0 in the 
SU(2) theory at L = 3 is shown in figurd?]. In this case the distribution is essentially zero 
away from a very sharp peak at cos a = 1 whose width is only 0(0.01). Notice, however, 
that this peak broadens as the continuum limit is approach and the expectation value of 
cos a subsequently falls. We have repeated the calculations of the supersymmetric Ward 

L < cos a > L < cos a > 

~3 0.86(8) ~3 0.95(2) 

j~ 0-86(3) ' ~4~ 0.79(5) 

~5~ 0.72(5) ~5~ 0.50(1) 

Table 14: SU(2) apbc fi = 0.25 Table 15: 517(2) apbc fi = 10.0 

identities for the thermal case and the results for QO\, Q12O4 and Q1O5 are shown in tables 
16, 17 and 18 at couplings // = 10.0 for SU(2) and fi = 5.0 for SU(3). For the thermal 
runs we are interested primarily in the data for large coupling /x which is now related 
to the inverse temperature - thus large [i corresponds to the theory at low temperature 
and/or large spatial volume. The data is also illustrated graphically in figure || for the 
case of SU(3). In this regime we find that the all the supersymmetric Ward identities 
computed within the phase quenched approximation are satisfied within statistical errors - 
thus we conclude that the use of thermal boundary conditions may indeed be quite useful 
in reducing phase fluctuation problems and eliminating the need for reweighting even away 
from the infinite volume limit. 

In the absence of spontaneous supersymmetry breaking expectation values computed 
in the thermal system should approach those of the periodic/zero temperature theory for 



A* 


Sb 


cR 
B 


0.25 


0.977(1) 


0.989(70) 


1.0 


0.981(1) 


0.995(45) 


5.0 


0.992(2) 


1.000(32) 


10.0 


0.999(2) 


1.003(160) 
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SU{2) 


SU{3) 


L 


B 


F 


B 


F 


3 


4.2(2) 


-4.2(2) 


5.2(3) 


-5.3(2) 


4 


1.8(2) 


-1.8(2) 


2.1(1) 


-2.1(1) 


5 


0.53(2) 


-0.55(2) 


0.90(8) 


-0.93(8) 



Table 16: QO x for SU(2) with apbc (jji = 10.0) and SU(3) (/i = 5.0) 





SU{2) 


SU(3) 


L 


B 


F 


B 


F 


3 


1.05(5) 


-1.03(6) 


1.30(6) 


-1.29(6) 


4 


0.45(5) 


-0.42(6) 


0.53(3) 


-0.50(4) 


5 


0.132(6) 


-0.101(6) 


0.22(2) 


-0.21(2) 



Table 17: Q 12 4 for SU{2) with apbc (// = 10.0) and SU{3) [p = 5.0) 





5C/(2) 


SU(3) 


L 




F 


B 


F 


3 


1.05(5) 


-1.03(6) 


1.30(6) 


-1.31(6) 


4 


0.45(5) 


-0.44(5) 


0.53(3) 


-0.51(3) 


5 


0.132(6) 


-0.119(5) 


0.22(2) 


-0.22(2) 



Table 18: Q1O5 for SU{2) with apbc (fj. = 10.0) and SU{3) (ji = 5.0) 

sufficiently low temperatures. No additional fine tuning should be required and this is seen 
in our data for the Ward identities. 

5.5 Rotational invariance 

To test for a restoration of rotational invariance as f3 is increased and the continuum limit 
is approached we have examined the two-point function 

G(x,y) =< B(0,0)B(x,y) > (5.7) 

where B(x,y) = [4>(x, y), 4>(x, y)\ and x and y are integer lattice coordinates. If the theory 
is rotationally invariant we would expect this correlator to be a function of just the radial 
distance r = \f x 2 + y 2 . Figure ||. shows a plot of this function for the SU (2) theory at 
(3 = 8.0 on a 8 x 8 lattice. The points lie within errors on a single curve lending support 
to the idea that at least the scalar sector of the theory is rotationally invariant at large 
distance. If all supersymmetries are restored for sufficiently small lattice spacing we would 
then expect the fermionic sector to also be rotationally invariant in the continuum limit. 
Fig pi] shows similar data for SU(3) at j3 = 4.0 and leads to similar conclusions. 
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5.6 Fluctuations of the scalar fields 

The classical vacua of this theory allow for any set of scalars which are constant over the 
lattice and satisfy 

[<M=0 (5.8) 

In the case of SU(2) and using the parameterization = fa + ifa we find vacuum states 
of the form 

01 = (A, 0,0) fa = (5,0,0) (5.9) 

together with global SU(2) rotations of this configuration. Thus we have a classical vacuum 
state for any value of A and B. In the case of SU(3) the general vacuum state can be 
realized by taking the fields fa and fa as constant diagonal matrices with arbitrarily large 
matrix elements. Generalization to arbitrary SU(N) the space of vacuum solutions is 
referred to as the classical moduli space. The presence of such a non-trivial moduli space 
corresponds to the existence of flat directions in the theory which have the potential to 
lead to divergences when inserted into a path integral. Notice, however, that the fermion 
operator develops zero modes precisely along these bosonic flat directions and thus quantum 
fluctuations of the fermions may suppress their contribution. In addition, entropic effects 
at large N may also play a role in removing potential divergences since the number of 
classical zero modes scales like N while the number of non-zero modes scales like N 2 . It 
then becomes a dynamical question as to whether the theory retains a non-trivial moduli 
space at the quantum level. 

To examine this issue we have measured the probability distribution of the scalar 



eigenvalues seen in the simulation. Figure 11. shows a plot of this distribution of scalar 
field eigenvalues for SU{2) obtained from high statistics runs (10 5 Monte Carlo trajectories) 
on a 3 x 3 lattice. We have mapped all data to positive values by exploiting the A — > —A 
symmetry of the action. The two sets of points correspond to (3 = 0.5 and (3 = 4.0. Contrary 
to the classical analysis it appears that the most probable scalar field configurations lie close 
to the origin in field space and this effect is enhanced as the coupling (3 is increased. A 
similar situation is seen for SU (3) in figure [12]. - the main difference being that there are 
now three peaks in this distribution (only two are seen after the mapping to the positive 
A axis) with a new peak appearing at the origin. This pattern repeats at larger TV - in the 
case of SU (N) we have observed that this probability distribution possesses N peaks |L7| . 
Since the scalars in this theory arise from dimensional reduction of gauge fields in four 
dimensional J\f = 1 Yang-Mills it seems likely that these N peaks in the scalar distribution 
are analogous to the N peaks appearing in the distribution of the expectation value of the 
Polyakov line in the usual deconfined phase of lattice gauge theory. 

At first glance these plots suggest that the scalars are driven to the origin as (3 — > oo 
(or equivalently L — > oo at fixed fj,) and the quantum continuum theory possesses only the 
trivial vacuum state < (f) >= 0. However, we must be be careful as the distributions clearly 
possess long tails extending out to large eigenvalue. We have examined this issue in more 
detail and indeed find good evidence for power law behavior in the tails of the distribution 



see figure 13. which shows a plot of log(P(A)) versus log A using data from the tail of 



the distribution for both SU(2) and SU(3) at a coupling (3 = 4.0. The exponent extracted 
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from a linear fit lies in the range of (minus) 2.1 — 2.3. Such an exponent would clearly lead 
to a normalizable distribution and hence a convergent partition function 4 . Moreover, since 
the global U(l) symmetry cannot break in two dimensions, it is clear that the first moment 
of this distribution < A > will vanish. However, a value of p < 3 would ensure that the 
variance < A 2 > — < A > 2 is divergent. This can be seen explicitly in our simulations 
with figure [14]. showing the Monte Carlo evolution of A 2 (averaged over the lattice and 
the number of colors) for the /3 = 4 SU(2) run. We see large fluctuations occurring at 
intervals of order 10, 000 RHMC trajectories. These spikes make it extremely difficult to 
determine the expectation value of < A 2 > - the statistical error does not decrease with 
the square root of the number of measurements as would be expected for a typical Monte 
Carlo process. These large spikes in < A 2 > are also seen for (3 = 0.5 and (3 = 8.0. 

We conjecture that the existence of these tails is related to the existence of a noncom- 
pact classical moduli space. This motivates us to consider the lattice theory in the limit of 
very large (3. In this limit we may restrict our attention to fields which are constant over 
the lattice and the theory reduces to a zero dimensional matrix model. This model has 
been studied previously with the result that the eigenvalue distributions again have power 
law tails with (negative) exponent p ~ 3 independent of the number of colors N 18]. 



The results for finite (3 in the two dimensional models look rather similar. It is tempting 
to interpret this divergence of the variance <A 2 > — <A> 2 as evidence that quantum 
effects will destroy any classical vacuum with < A >^ - the fluctuations around any 
such solution will be so large as to swamp the mean value and lead to a restoration of the 
A — > — A symmetry. 

6. Conclusions 

In this paper we have extended our previous numerical study of N = 2 super Yang-Mills 



theory in two dimensions described in [13]. The discretization we have employed was first 
proposed in Q and is based upon mapping the twisted continuum theory into a theory 
of Kahler-Dirac fields which may then be discretized preserving gauge invariance and a 
single supersymmetry while simultaneously avoiding problems of fermion doubling. In the 
previous paper we presented numerical results supporting the existence of an exact scalar 
supersymmetry Q in the case of the SU{2) theory. In this paper we have derived the 
action of the remaining twisted supersymmetries on the component Kahler-Dirac fields 
and studied some non-trivial Ward identities following from those symmetries for both 
SU(2) and SU(3). It should be noted that many of these additional Ward identities are 
automatically zero because of gauge invariance and so the additional symmetries supply 
rather few additional constraints. 

Our numerical results, carried out within the phase quenched approximation, support 
a restoration of all supersymmetries in the continuum limit at least for large continuum 
dimensionless coupling \x = g 2 A. Such a limit corresponds to large physical volumes and 
appears necessary in order to justify the phase quenched approximation which we have, of 



4 Actually since here we use periodic boundary conditions on all fields the partition function is more 
properly thought of as a Witten index. 
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necessity, employed. Furthermore, we have shown that an additional suppression of this 
phase on finite lattices can be achieved by using thermal boundary conditions. We have 
also examined a correlation function of the scalar fields which shows good evidence of a 
restoration of rotational invariance at large distances which is also encouraging. 

These conclusions are consistent with perturbative arguments based on continuum 
power counting see for example [10, which argue that a single exact supersymmetry 
plus gauge invariance and the additional global U(l) symmetry prohibit the occurrence of 
relevant SUSY violating operators as the cut-off is removed. 

Finally, we present new results concerning the eigenvalue distributions of the scalar 
fields. For SU (2) we see two symmetrically disposed peaks in this distribution which 
narrow and approach the origin as (5 — > oo. In the case of SU{3) three peaks are seen. 
This pattern appears to persist at larger N with the distribution showing N peaks. A more 
detailed analysis reveals power law tails to these distributions which ensures that sufficiently 
high moments of the eigenvalue distributions will diverge. Our data is consistent with a 
divergence of the second moment for both SU(2) and SU(3) at least for large (5. Such a 
result is similar to previous results for supersymmetric matrix models [18]. 

Furthermore, we have argued that the divergent variance would lead to large fluctu- 
ations in the scalar fields which would tend to wash out any classical vacuum state with 
< 4> 0. Thus the classical moduli space does not survive in the full quantum theory. 
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Figure 2: QO x vs L for SU(3) with fi = 0.25,5.0 
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SU(3) 




Figure 3: Q 12 4 vs L for SU(3) with fj, = 0.25,5.0 
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cos(alpha) 

Figure 5: P(cosa) vs cos a for SU{2) with L = 3 fi = 0.25 
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SU(2) mu=10.0 




cos(alpha) 

Figure 6: P(cosa) vs cos a for SU{2) with L = 3 ^ = 10.0 
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SU(2) mu=10.0 
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Figure 7: P(cosa) vs cos a for thermal SU(2) with L = 3 [i = 10.0 
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Thermal SU(3) 





Figure 9: G(r) vs r for SU(2) and = 8.0 with L = 8 
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Figure 10: G{r) vs r for SU(%) and (3 = 8.0 with L = 
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Probability distribution of eigenvalues of scalars 
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Figure 11: Eigenvalue distribution for SU(2) and L = 3 
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Probability distribution of eigenvalues of scalars 

SU(3),3 2 
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Figure 12: Eigenvalue distribution for ST/ (3) and L = 3 
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Figure 13: log (P(A) vs log A for L = 3 and both SU{2) and SU(3) 
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SU(2), beta=4.0 
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Figure 14: Monte Carlo history for L = 3, (3 = 4.0 and SU{2) 
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